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Abstract 

Null results of experimental searches for the Higgs boson and the superpartners imply a 
certain amount of fine-tuning in the electroweak sector of the Minimal Supersymmetric 
Standard Model (MSSM). The "golden region" in the MSSM parameter space is the 
region where the experimental constraints are satisfied and the amount of fine-tuning 
is minimized. In this region, the stop trilinear soft term At is large, leading to a 
significant mass splitting between the two stop mass eigenstates. As a result, the 
decay I2 — ► t\Z is kinematically allowed throughout the golden region. We propose 
that the experiments at the Large Hadron Collider (LHC) can search for this decay 
through an inclusive signature, Z + 2% + Ifsp + X. We evaluate the Standard Model 
backgrounds for this channel, and identify a set of cuts that would allow detection of 
the supersymmetric contribution at the LHC for the MSSM parameters typical of the 
golden region. We also discuss other possible interpretations of a signal for new physics 
in the Z + 2% + frp + X channel, and suggest further measurements that could be used 
to distinguish among these interpretations. 



1 Introduction 



It is widely believed that physics at the TeV scale is supersymmetric. The simplest realistic 
implementation of this idea, the minimal supersymmetric standard model (MSSM) [1, 2], 
is the most popular extension of the standard model (SM). However, null results of ex- 
perimental searches for the superpartners and, especially, the Higgs boson, place non-trivial 
constraints on the parameters of the model. Furthermore, the requirement that the observed 
electroweak symmetry breaking (EWSB) occur without significant fine-tuning places an ad- 
ditional constraint. It is well known that there is a certain amount of tension between these 
two constraints [3]. Several authors have interpreted this tension as a motivation to extend 
the minimal model [4], or to question the conventional ideas about naturalness [5]. An alter- 
native interpretation, which we will explore in this paper, is that data and naturalness point 
to a particular "golden" region within the parameter space of the minimal model, where the 
experimental bounds are satisfied and fine-tuning is close to the minimum value possible 
in the MSSM. This minimal value itself depends on the messenger scale of supersymmetry 
breaking A mess , determined by dynamics outside of the MSSM, in addition to the MSSM 
parameters 1 . However, for any A mess , the points in the golden region require less fine-tuning 
compared to the rest of the MSSM parameter space. Thus, independently of the model of 
SUSY breaking, nature seems to provide us with a hint about what the MSSM parameters 
might be 2 . In this paper, we will discuss experiments at the Large Hadron Collider (LHC) 
which will be able to determine whether this hint is correct. 

Both the Higgs mass bound and naturalness considerations probe the effective Higgs 
potential, which is primarily determined by the parameters of the Higgs and top sectors of 
the MSSM. (The quantum part of the potential is dominated by the top/stop loops due to a 
large value of the top Yukawa coupling.) It is therefore these sectors that are most directly 
constrained by data. We will focus on collider measurements probing these sectors. 3 

The golden region is characterized by relatively small values of the fi parameter and the 
stop soft masses mg3, m u 3 (both are required to minimize fine-tuning of the Z mass), and 
a large stop trilinear soft term A t (required to raise the Higgs mass above the LEP2 lower 
bound). The spectrum is then expected to contain light neutralinos and charginos with a 
substantial higgsino content, as well as two light (sub- TeV) stop mass eigenstates, t± and £2, 
with a large (typically a few hundered GeV) mass splitting. A striking consequence of such 
a "split stop" spectrum is that the decay 

t 2 -> fi + Z (1) 

1 For example, it was claimed in Ref. [6] that in models with "mirage mediation" of SUSY breaking [7] 
the scale A mess can be as low as 1 TeV, resulting in fine-tuning of 20% or better. See Rcf. [8] for a discussion 
of difficulties in realizing such a scenario, and Ref. [9] for an alternative implementation. 

2 An explicit model of supersymmetry breaking in a grand-unified framework which naturally generates 
SUSY breaking parameters in the golden region was constructed in [10]. 

3 Our approach is more model-independent than that of Kitano and Nomura in Ref. [11], which also consid- 
ered collider signatures of the MSSM with parameters in the golden region. However, the signatures studied 
in [11] mainly probe the features of the superpartner spectrum dictated by a specific (mirage mediation) 
model of SUSY breaking [6], rather than the direct consequences of data and naturalness. 
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is kinematically allowed. Observing this decay at the LHC would provide clear evidence 
that the stop mass difference is larger than the Z mass, and studying the Z distributions 
would provide an approximate measurement of this quantity. In this paper, we will argue 
that the decay (1) should be observable at the LHC, with realistic integrated luminosity, for 
the MSSM parameters in the golden region. 

The experimental signature of the decay (1) depends on the decay pattern of the t\. Since 
stops are almost always pair-produced at the LHC, it also depends on how the second ti 
decays. The details of both decay patterns depend on the superpartner spectrum. However, 
both t\ and ti decay products always contain a b quark, produced either directly or through 
a top decay, as well as (under the usual assumptions of conserved R parity and weakly 
interacting lightest supersymmetric particle) large missing transverse energy. We therefore 
propose an inclusive final state 

Z + 2j b + fir + X, (2) 

where Z is assumed to be reconstructed from leptonic decays and j b denotes a b jet, as a 
signature of the £2^2 production followed by the decay (1). 

Throughout the golden region of the MSSM, both the I2 pair-production cross section 
and the branching fraction of the decay (1) are sizeable. Therefore, a null result of a search 
for a non-SM contribution in the channel (2) would provide a strong argument against this 
scenario. Unfortunately, a positive identification of non-SM physics in this channel would 
not necessarily imply that the stops are split. Indeed, in the MSSM, events in this channel 
may appear even if the decay (1) is kinematically forbidden, since Z bosons may also be 
produced in decays of neutralinos and charginos [12]. For example, a cascade 

b^bxl x°->Zx?, (3) 

or a similar cascade with charginos replacing the neutralinos, gives the signature (2). Dis- 
tinguishing these interpretations is difficult, and there is no single "silver bullet" observable 
that would remove this ambiguity. However, a variety of measurements can be used to shed 
light on this question (see Section 5), and combining all available evidence may allow one 
to build a convincing case for (or against) the interpretation of the signature (2) in terms of 
the decay (1). 

The paper is organized as follows. In Section 2, we review the fine-tuning and Higgs mass 
constraints in the MSSM, as well as other experimental results that determine the shape of 
the golden region. In Section 3, we define a benchmark point which is characteristic of the 
golden region and suitable for studying its collider phenomenology. Section 4 is dedicated to a 
detailed analysis of the observability of the Z+2j b +fir signature, including a study of the SM 
backgrounds. In Section 5, we discuss the alternative interpretations of this signature within 
the MSSM, and outline the measurements that would need to be performed to discriminate 
between these interpretations. Section 6 contains our conclusions, and outlines some possible 
directions for future work. 
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2 The Golden Region 



In this section, we will discuss the constraints on the MSSM parameters imposed by current 
experimental data and naturalness, focusing on the Higgs and top sectors. Our goal is to 
understand the qualitative features of the MSSM golden region, rather than to determine 
the precise location of its boundaries which are in any case fuzzy due to an inherent lack of 
precision surrounding the concept of fine-tuning. With this motivation, we will make several 
approximations which greatly clarify the picture. 

Phenomenological studies of the MSSM are complicated by the large number of free 
parameters. Typically, studies are performed within simplified frameworks, which assume 
certain correlations among the parameters motivated by high-scale unification and/or by 
specific models of SUSY breaking. However, the shape of the golden region is to a great 
extent independent of such assumptions. The Higgs sector of the MSSM is strongly coupled 
to the top sector, but couplings to the rest of the MSSM are weaker. One may therefore 
begin by considering the Higgs and top sectors in isolation; that is, the gauge and non-top 
Yukawa couplings are set to zero. In this approximation, physics is described in terms of the 
holomorphic Higgs mass /i and the six parameters appearing in the soft Lagrangian for the 
Higgs and top sectors: 

C = -m2|#J 2 -m2|#,| 2 -(M^^ 

(4) 

where yt is the MSSM top Yukawa coupling, y t = yf M /sinj3. Since the model has to 
reproduce the known EWSB scale, v = 174 GeV, only six parameters are independent. We 
choose the physical basis: 

t&nP,[i,m A ,m 1 ,m 2 ,9 t , (5) 

where uia is the CP-odd Higgs mass, rh\ and rh 2 are stop eigenmasses (by convention, 
m 2 > mi) and 6 t is the stop mixing angle. We will analyze the fine-tuning and Higgs 
mass constraints in this approximation and map out the golden region in the six-parameter 
space (5). 

Before proceeding, let us discuss the sizes of contributions to the relevant observables that 
are omitted in this approximation scheme. The leading contributions to the Higgs effective 
potential due to the SU(3), SU(2) and U(l)y gauge interactions and the bottom Yukawa 
coupling 4 are expected to be of the order 

gjMi g\Ml gjMl m\Ml tan 2 (3 
167r 2 A 2 Mf \\Mf A 2 Mf m 2 M| ' K) 

respectively, compared to the one-loop top sector contribution. Here gi and Mj are the gauge 
couplings and (weak-scale) gaugino masses for each group, and Mi is the stop mass scale, 
which can be conveniently taken as the average between the two stop eigenmasses. (The 

4 The corrections due to other Yukawa couplings are always negligible. 
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Figure 1: Contours of 1% fine-tuning in the (/i, tan/3) plane. The black (solid) contour 
corresponds to = 100 GeV, but remains essentially unchanged for any value of in the 
range between 100 and 1000 GeV. The red (dashed) and blue (dotted) contours correspond 
to m A = 1.5 and 2 TeV, respectively. 



same definition can be made for M~ b if sbottoms are non-degenerate.) For a wide range of 
sensible superpartner spectra, these corrections are subdominant: this is the case if 



M 1 /M i <4, M 2 /Mi<2, M 3 /Mt< 10, M- b < 



35Mt 
tan (3 



(7) 



The following discussion is valid for spectra obeying these constraints. If some of the above 
inequalities are violated, the analysis could be easily extended to include the corresponding 
effects; however, little additional insight would be gained. 



2.1 Constraints on the Higgs Sector 

At tree level, the Z mass in the MSSM is given by 
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Following Barbieri and Guidice [13], we quantify fine-tuning by computing 
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where £ = m 2 u ,m\,b, \i are the relevant Lagrangian parameters. In terms of the physical 
parameters (5), we obtain 
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where we assumed tan/5 > 1. The overall fine-tuning A is defined by adding the four A's 
in quadruture; values of A far above one indicate fine-tuning. For concreteness, we will 
require A < 100, corresponding to fine tuning of 1% or better. This requirement maps out 
the golden region in the space of (tan/3, Ma), as illustrated in figure 1. (We do not plot 
/i < 100 GeV, since this region is ruled out by LEP2 chargino searches.) The shape of this 
region is easily understood. In the limit of large tan/5, the parameters A(m 2 ) and A(m^) 
are small, and A(/i) and A(b) (considered separately) lead to constraints 



m z 



A 1 / 2 m A A 1 / 2 
2 mz 2 



(12) 



which are clearly reflected in Fig. 1. As f3 approaches 7r/4, the factors of 1/ cos 2/3 and 
tan 2/3, present in all four A parameters, become large, and as a result the model is always 
fine-tuned for tan/5 ^ 2. 



2.2 Constraints on the Top Sector 

Naturalness also constrains the size of the quantum corrections to the parameters in Eq. (8). 
The largest correction in the MSSM is the one-loop contribution to the m 2 parameter from 
top and stop loops: 

3 2A 2 
6m Hu ~ (^Qs + ^4) +^ 2 (Asin/5-/icos/5) 2 ) log 

« (rf (*! + *3 - 2m?) + ^M>! sin * 2S( ) log _^!_ (13) 

where m t is the top mass, A is the scale at which the logarithmic divergence is cut off, and 
finite (matching) corrections have been ignored. The correction induced by this effect in the 
Z mass is 

6,ml « -*< (l - - J^) , (14) 
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where we ignored the renormalization of the angle (3 by top/stop loops: the contribution of 
this effect scales as 1/ tan 2 /? and is subdominant for tan/5 ^ 2. To measure the fine-tuning 
between the bare (tree-level) and one-loop contributions, we introduce 



8 t m% 



(15) 



Choosing the maximum allowed value of A t selects a region in the stop sector parameter 
space, (rhi,m2,0 t ), whose shape is approximately independent of the other parameters. 5 
This constraint is shown by the black (dashed) lines in Figs. 2, where we plot 5%, 3%, 1% 
and 0.5% tuning contours (corresponding to A t = 20, 33.3, 100, and 200, respectively) in the 
stop mass plane for several values of 9 t and tan (3 = 10. Note that the particular values of 
A t depend on the scale A; we choose it to be 100 TeV in this figure. However, the shape 
of the contours and the obvious trend for tuning to increase with the two stop masses is 
independent of A. 

The second constraint that determines the shape of the golden region is the LEP2 lower 
bound on the Higgs mass [14]. For generic MSSM parameter values, the limit on the lightest 
CP-even Higgs is very close to that for the SM Higgs: 

m(h°) > 114 GeV. (16) 

It is possible for a lighter Higgs (down to about 90 GeV) to be consistent with the negative 
results of the LEP2 searches; however, this requires precise coincidence between m(h°) and 
rriA, which should be regarded as additional source of fine-tuning. Thus, we will use the LEP2 
bound for the SM Higgs [15], 114.4 GeV, as the lower bound on m(h°) in this analysis. At 
tree level, the MSSM predicts m(h°) < mz \ cos2/5|, and large loop corrections are required 
to satisfy this bound. Extensive calculations of these corrections have been performed in the 
literature (for a recent summary of the status of these calculations, see Ref. [16]). Complete 
one-loop corrections within the MSSM are known. The dominant one-loop contribution is 
from top and stop loops; for tan (3 ^ 35, the sbottom loop contribution is also important. 
The two-loop corrections to these contributions from strong and Yukawa interactions are 
also known. Numerical packages incorporating these results are available [17, 18]. For our 
purposes here, however, it is convenient to use a simple analytic approximation, due to 
Carena et. al. [19], which includes the one-loop and leading-log two-loop contributions from 
top and stop loops: 

2 . 



m 2 (/i°) = m 2 z cos 2 2(3 1 
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5 Note that we choose not to combine the tree- level and quantum fine-tuning measures into a single tuning 
parameter; doing so would make the analysis less transparent without producing additional physical insights. 
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Figure 2: Fine-tuning (black/dashed contours), Higgs mass bound (red/colid contours), and 
p-parameter (blue/dotted contours) constraints in the (mi, 5m) plane. The six panels corre- 
spond to (starting from the upper-left corner, clockwise): 9 t = 0, 7r/25, vr/15, tt/6, tc/ 4, it/3. 
In all panels tan/3 = 10. The yellow/shaded intersection of the regions allowed by the three 
constraints is the MSSM "golden" region. 



where a 3 is the strong coupling constant evaluated at the pole top quark mass M t ; m t 
M t /(1 + ^^3) is the on-shell top mass; and 

2(A- / ucot/3) 2 ( {A t - fi cot /3) 2 \ 

stt. — - 



M 2 \ 12M 2 / ' 

susy \ susy / 

M 2 

t = log^. (18) 

The scale M s 2 usy is defined as the arithmetical average of the diagonal elements of the stop 
mass matrix. The expression (17) is valid when the masses of all superparticles, as well as 
the CP-odd Higgs mass m^, are of order M susy . Additional threshold corrections may be 
required, for example, if rriA < M susy ; for simplicity, we will ignore such corrections here. 
Eq. (17) agrees with the state-of-the-art calculations to within a few GeV for typical MSSM 
parameters [16]; while such accuracy is clearly inadequate for precision studies, it is sufficient 
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for the present analysis. 6 

The contours in the stop mass plane corresponding to the LEP2 Higgs mass bound 
are superimposed on the fine-tuning contours in Figs. 2. The positions of these contours 
depend strongly on the top quark mass. We used M t = 171.4 ±2.1 GeV [20], and plotted the 
constraint corresponding to the central value (thick red/solid lines), as well as the boundaries 
of the 95% c.l. band (thinner red/solid lines). The contours are approximately independent 
of tan/5 for 3 ^ tan/3 ^ 35; the golden region shrinks rapidly outside of this range of tan/3. 
We use tan f3 — 10 in the plots. The overlap between the regions of acceptably low fine- 
tuning (for definiteness, we choose A t = 100) and experimentally allowed Higgs mass defines 
the golden region, shaded in yellow in Figs. 2. 

2.3 Collider Bounds, Precision Electroweak Constraints and Rare 
Decays 

Apart from the Higgs mass bound, several other observables constrain the shape of the 
golden region. 

First, direct collider bounds play a role in determining the boundary at low p and fh\. 
LEP2 searches for direct production of charginos and stops constrain both p and fh\ to be 
above ~ 100 GeV, and are to a large extent independent of the rest of the MSSM parameters. 
(At large tan/3, it can be easily shown that m(xf) < \p\ for any M 2 .) The Tevatron stop 
searches yield a similar (though more model-dependent) bound on mi. A sbottom search in 
the bxi channel (which is relevant because the b\ mass is given by mg3, and can be expressed 
in terms of fhi and m 2 ) places a lower bound m(b L ) > 200 GeV. However, this bound will 
not be used in our analysis since it is highly sensitive to the neutralino mass and can be 
easily evaded if m(xi) > 80 GeV. 

Second, in the presence of a large A t term, stop and sbottom loops may induce a signifi- 
cant correction to the p parameter. This correction is known at the two- loop level [21]; for 
our purposes, it suffices to use the one-loop result: 

3G 

Ap = (- sin 9 2 t cos 6*F (ml, m 2 ) + cos 2 ^F (m 2 , m\ L ) + sin 2 6 t F (m 2 2 , mfj) , (19) 

where 

FoM) = a + 6-^-log^. (20) 
a — b b 

Expressing rag in terms of mi, m 2 and 9 t , and using the PDG value p = 1.0002 io;ooo7 [15], 
we obtain the 95% c.l. contours in the stop mass plane shown by the blue/dotted lines in 
Figs. 2. This constraint eliminates a part of the parameter space with very low mi and large 
5m. 

6 Wc also verified that the Higgs mass at the benchmark point used for the collider phenomenology analysis 
in this paper satisfies the LEP2 bound with a more precise numerical calculation using SuSpect; see Section 
3. 
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Table 1: The benchmark point: MSSM input parameters, defined at the weak scale. (All 
dimensionful parameters are in GeV.) 



Finally, several low-energy measurements play a role in constraining the MSSM parameter 
space; among these, the b — > decay rate [22] and the anomalous magnetic moment of the 
muon, — 2 [23] , provide the most stringent constraints. The supersymmetric contribution 
to — 2 depends sensitively on the slepton and weak gaugino mass scales, and only weakly 
on the parameters defining the golden region. On the other hand, since the golden spectrum 
contains light stops and higgsinos, we can expect a large contribution to the b — > 37 rate 
from the t — H loop. It is well known, however, that this can be cancelled by the contribution 
of the top-charged Higgs loop. A simplified analysis of this constraint based on the one-loop 
analytic formulas presented in Ref. [24] shows that for any values of the stop masses inside 
the golden region in Figs. 2, and for any value of /i between 100 and 500 GeV, one can find 
values of in the 100-1000 GeV range for which this cancellation ensures consistency with 
experiment. (Recall that m 2 (iJ ± ) = m^ + m^.) For low rh\ and /i, however, the cancellation 
only occurs in a narrow band of tjia, which can be thought of as an additional source of fine 
tuning. A detailed analysis of this issue is outside the scope of this paper. 



3 A Benchmark Point for Collider Studies 

The analysis of Section 2 defined the golden region in the six-dimensional parameter space (5); 
its shape is approximately independent of the other MSSM parameters. This region has the 
following interesting qualitative features: 

• Both stops typically have masses below 1 TeV; 

• A substantial mass splitting between the two stop quarks is required: typically, 
5m > 200 GeV; 

• The stop mixing angle must be non-zero: there is no intersection between the natural- 
ness and Higgs mass constraints for 9 t = 0, ir/2. 

The first feature implies that both ti and t 2 will be produced with sizeable cross sections at 
the LHC, so that the stop sector can be studied directly experimentally. The second feature 
implies that the decay mode t 2 — > t±Z is kinematically allowed. The vertex responsible for 
this decay is given by 



2 Sm2ft tl2-3<)' (21 > 
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Table 2: The benchmark point: physical spectrum. All masses are in GeV. The masses of 
the superparticles not listed here are close to 1 TeV. 
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Table 3: The benchmark point: branching ratios of t 2 decay modes, in %. 



where c w and s w are the cosine and sine of the SM Weinberg angle. The last of the three 
points above then guarantees that the vertex is non-zero, and the decay t 2 — > t\Z indeed 
occurs. 

The branching ratio of the t 2 —> t\Z mode depends on which competing t 2 decay channels 
are available. The possible two-body channels are 

tg, fx , bx+, bW + , bH + , ^h , hH , ^A , (22) 

where x° an d X + denote all the neutralinos and charginos that are kinematically accessible, 
and flavor- changing couplings are assumed to be negligible. 

We would like to evaluate the prospects for observing the t 2 — > l\Z decay mode at the 
LHC. For concreteness, we choose a benchmark point (BP) within the golden region, and 
perform a detailed analysis of the signal at this point (see Section 4). The BP is defined 
in terms of the weak-scale MSSM parameters. We assume that all soft parameters are 
flavor-diagonal. Further, we assume a common soft mass for the first and second generation 
squarks, rriq = mq\,i = m u i,2 = m d i,2, and for all sleptons, = m L i,2,z = m e i,2,3 = 77^1,2,3. 
All A terms have been set to zero, with the exception of A t . The parameters defining the 
BP are listed in Table 1. The top and Higgs sector parameters are chosen so that the BP is 
comfortably inside the golden region, well away from the boundaries, and is representative 
of this region. In particular, the lightest Higgs mass at the BP is well above the LEP bound. 
(The physical spectrum of the model at the BP was computed using the SuSpect software 
package [18] and is listed in Table 2.) Gaugino, slepton, and first and second generation 
squark masses are set at 1 TeV. Varying these parameters does not have a significant effect 
on the stop production rate and decay patterns, and thus the conclusions of the analysis 
in Section 4 are largely independent of these choices. Using SuSpect, we checked that 
the b — > S7 branching ratio, the p parameter and the supersymmetric contribution to the 
muon anomalous magnetic moment at the BP are consistent with the current experimental 
constraints. 
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Figure 3: Cross section of the process pp — > it* at the LHC, — 14 TeV, at tree level. 
Factorization and renormalization scales were set to /i — Mj, and the CTEQ6L1 parton 
distribution function set [27] was used. 



The t 2 decay branching ratios at the BP were evaluated using the SDECAY package [25], 
and are listed in Table 3. The t 2 - ► t\Z mode has a substantial branching ratio, about 
31%. Note that of the possible t 2 decay modes listed in Eq. (22), only the tg channel is 
kinematically forbidden at the BP. If the gluino mass were lowered to allow this decay, 
the branching ratio of the t 2 — > t\Z mode would be suppressed. However, this effect is 
not dramatic: we checked that if M3 is varied between 300 and 1000 GeV, keeping all other 
MSSM parameters fixed at their values listed in Table 1, we still obtain Br(t 2 — > i\Z) 17%. 
This is an example of the robustness of the stop decay pattern with respect to the variations 
of the non-stop sector MSSM parameters, mentioned above. 



4 Observability of the Z + 2jb + fk + X Signature at the 



Stop pair production cross section at the LHC, computed using the MadGraph/MadEvent v4 
software package [26], is shown in Fig. 3. At the benchmark point, we find a(pp — > t 2 t 2 ) = 
0.05 pb, corresponding to about 500 t 2 pairs per year at the initial design luminosity of 10 
fb~ 1 /y ear - The produced t 2 decays promptly, with branching ratios listed in Table 3; in 
about 52% of the events, either one or both of the produced stops decays in the t x Z mode. 
This decay is followed by a cascade 
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(23) 
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where the jets and leptons produced in the xt decays are very soft due to a small chargino- 
neutralino mass splitting. The details of this cascade are particular to the chosen BP, and are 
quite model-dependent. There are, however, two model-independent features true for all ti 
and t\ decays: the cascade always contains a b jet (produced either directly or via top decay), 
and it always ends with the lightest supersymmetric particle (LSP), the neutralino Xi, giving 
a missing transverse energy signature. In order to make the analysis as model-independent 
as possible, we focus on an inclusive signature, 

Z(£ + ,£-) + 2j b + fa + X, (24) 

where Z(£ + ,£~) denotes a lepton pair (£ = e or fi) with the invariant mass at the Z peak. 
The presence of energetic leptons ensures that essentially all such events will be triggered on; 
we will assume a triggering probability of 1 for this analysis. Note that events with hadronic 
Z decays may in principle be triggered on due to large fa; however, this sample would 
suffer from a severe background of purely QCD events with apparent fa due to jet energy 
mismeasurement, and it will not be used in our study. Note also that the requirement that 
both jets be 6-tagged can be relaxed, as will be discussed below. While a cleaner sample is 
obtained if two 6-tags are required, this sample is smaller due to the less-than-perfect tagging 
efficiency, which may be relevant since the signal rates are not large. 

To assess the observability of the signature (24), we have simulated a statistically signifi- 
cant event sample for the signal and several SM background channels 7 using the MadGraph/ 
MadEvent v4 software package [26]. This tool package allows us to generate both SM and 
MSSM processes, so that the signal and backgrounds can be treated uniformly. The parton 
level events generated by MadEvent were recorded in the format consistent with the Les 
Houches accord [28, 29]. These events were then passed on to the Pythia package [30], 
which was used to simulate showering and hadronization, as well as the decays of unstable 
particles. Finally, the Pythia output was processed by the PGS 3.9 package [31], which 
provides a simple and realistic simulation of the response of a "typical" particle detector. 
(A more detailed analysis of the detector effects using complete ATLAS and CMS detector 
simulation packages would clearly be interesting, but is outside the scope of this study.) 
The final output was analyzed with ROOT, using only detector level information for event 
reconstruction. 

The following SM backgrounds have been considered in detail: 

• jjZZ , which can produce the signature (24) if one Z decays invisibly and the other 
one is reconstructed in £ + £~; 

• tiZ, with Z — > £ + £~ and one or both tops decaying leptonically (with fa due to neutri- 
nos), or both tops decaying hadronically (with fa due to jet energy mismeasurement). 

7 At the chosen benchmark point, the events containing i 2 — * i\Z arc the only non-SM source of the 
signature (24), so there arc no "signal backgrounds". Possible alternative interpretations of this signature 
in the general MSSM context are discussed below in Section 5. 
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signal: t2^2 


33ZZ 


4-4- V 

ttz 


tt 


A. A. V 


Oprod(pb) 


0.051 


0.888 


0.616 


552 


824 


total simulated 


9964 


159672 


119395 


3745930 


1397940 


1. leptonic Z(s) 


1.4 


4.5 


2.6 


0.04 


2.1 


2(a). Pt(ji) > 125 GeV 


89 


67 


55 


21 


41 


2(b). p t (j 2 ) > 50 GeV 


94 


93 


92 


76 


84 


3. o-tag 


H A 

64 


8 


A A 

44 


57 


5 


4. 7 (Z) > 2.0 


89 


66 


69 


26 


68 


5. pr > 225 GeV 


48 


2.2 


4.4 


1.7 


< 0.9 (95% c.l.) 
(ext.) 


iV exp (100 fb- 1 ) 


16.4 


2.8 


10.8 


8.8 


< 177 (95% c.l.) 

(ext.) 



Table 4: Summary of the analysis of observability of the supersymmetric golden region 
signature (24). First row: Production cross section for the signal and background processes 
at the LHC. Second row: Number of Monte Carlo events used in the analysis. Rows 3-8: 
Cut efficiencies, in%. Last row: The expected number of events for an integrated luminosity 
of 100 fb" 1 . 



• tt, with both tops decaying leptonically and the invariant mass of the two leptons 
accidentally close to mz- 

The total production cross sections (with p™^ t = 50 GeV for the jjZZ channel) and the 
size of the event sample used in our analysis for each channel are listed in the first two rows 
of Table 4. To identify the events matching the signature (24), we impose the following set 
of requirements on the event sample: 

1. Two opposite-charge same-flavor leptons must be present with ^Js(£ + £~) = M z ± 2 
GeV. 

2. Two hard jets must be present, with pt > 125 GeV for the first jet and px > 50 GeV 
for the second jet; 

3. At least one of the two highest-pr jets must be 6-tagged; 

4. The boost factor of the Z boson, ^y(Z) = l/\Jl — f|, reconstructed from the lepton 
pair, must be larger than 2.0; 

5. A missing E T cut, $t > 225 GeV. 

The efficiencies of these cuts are given in Table 4, and the pr distribution of the events 
passing cuts 1-4 is shown in Figure 4. While the overall rate of the SM background processes 
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Figure 4: Missing distribution of the events passing cuts 1-4. Signal is shown in black; 
jjZZ, tiZ and ti backgrounds are shown in blue/dark-gray, green/gray, and yellow/light- 
gray, respectively. The normalization corresponds to an integrated luminosity of 300 fb^ 1 at 
the LHC. 



is much higher than the signal rate, the cuts 1-5 are quite effective in discriminating signal 
from background. Assuming that the search is statistics-limited, we estimate that a 3-sigma 
observation would require 75 fb _1 of data, while a definitive 5-sigma discovery is possible 
with 210 ft) -1 . Note that one important contribution to the background, from the tt channel, 
can be effectively measured from data by measuring the event rates with dilepton invariant 
masses away from the Z peak and performing shoulder subtraction. This procedure is likely 
to be statistics-limited. However, systematic uncertainties in other background contributions 
could play a role in limiting the reach, and should be studied carefully with a more detailed 
detector simulation. 

We also briefly considered several other irreducible SM backgrounds which are expected 
to be less significant than the ones listed in Table 4, but might nevertheless be relevant. 
The most important one of these is tij, where j is a hard jet. The cross section for this 
channel is suppressed compared to ti, but the presence of the additional hard jet increases 
the probability that the events will pass the jet pt cut (cut 2). We find a parton- level cross 
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section a(tij,p? T > 125 GeV) = 65 pb. Assuming conservatively that all these events pass 
the cut 2, and that the efficiencies of all other cuts are the same as for the tt sample, we 
expect that this background would add at most about 50% to the tt rate. As in the tt case, 
this contribution can be subtracted using data away from the Z peak in the lepton invariant 
mass distribution. Assuming that the statistical error dominates this subtraction, the net 
effect would be an increase in the integrated luminosity required to achieve the same level 
of significance by at most about 10%. 8 Other backgrounds we considered are three vector 
boson channels ZZZ, ZZW , and ZWW; as well as channels with single top production, tZj 
and tZj. Combining the parton- level cross sections for these channels with the branching 
ratios of decays producing the signature (24) results in event rates that are too small to 
affect the search. 

While the SM processes considered above genuinely produce the signature (24), other 
SM processes may contribute to the background due to detector imperfections. We expect 
that the dominant among these is the process jjZ, with Z — > and apparent fir due 
to jet energy mismeasurement or other instrumental issues. We conducted a preliminary 
investigation of this background by generating and analyzing a sample of 1.4 x 10 6 jjZ 
events with p T u ^, t = 50 GeV (see the last column of Table 4). None of the events in this 
sample pass the cuts 1-5. This allows us to put a 95% c.l. bound on the combined efficiency 
of this set of cuts for the jjZ sample of about 2 x 10~ 6 , corresponding to a background 
rate about 10 times larger than the signal rate. However, we expect that the actual jjZ 
background rate is well below this bound, since all 349 events in our sample that pass the 
cuts 1-4 in fact have fir below 50 GeV. We find that the fir distribution of these 349 events 

can be fit with an exponential, N oc e~ ai0 -^ r ', where fip is in units of GeV. Assuming that 
this scaling adequately describes the tail of the distribution at large fir, we estimate that the 
rate oijjZ events passing all 5 cuts is completely negligible and that this background should 
not present a problem. This conclusion is of course rather preliminary, and this issue should 
be revisited once the performance of the LHC detectors is understood using real data. Note 
that the necessity to understand the shape and normalization of the large apparent fir tail 
from SM processes with large cross sections is not unique to the signature discussed here, 
but is in fact crucial for most SUSY searches at the LHC. 

As an alternative, we considered a variation of the analysis where the cuts 1, 2, and 5 
are unchanged, cut 4 is eliminated, and two 6-tagged jets are required. The cut efficiencies 
for this analysis are summarized in Table 5. Unfortunately, the Monte Carlo samples used 
in our analysis are not large enough to reliably estimate the efficiencies of this set of cuts 
applied to the backgrounds, since only one event out of all background samples passes the 
cuts. Therefore we list the 95% c.l. upper bounds on the efficiencies, and on the number of 
background events expected for a 100 ftr 1 event sample, in the table. It is clear that while 
the second 6-tag is quite efficient in improving the S/B ratio, this search suffers from low 

8 The ttj background may be suppressed very effectively by requiring that the hardest jet be 6-tagged (as 
oppossed to one of the two hardest jets in our main analysis), since the extra jet is always initiated by a 
gluon or a light quark. However this would also reduce the signal and all other backgrounds by about a half 
due to the lower probability of tagging a single jet, resulting in lower significance. 



15 





signal: £2^2 




4-4- V 

ttz 


tt 


A. A. V 


Oprod(pb) 


0.051 


0.888 


0.616 


552 


824 


j j "I ■ 1 j 1 

total simulated 


9964 


159672 


119395 


3745930 


1397940 


1. leptonic Z(s) 


1 A 

1.4 


A £ 

4.0 


p. 
z.o 


n n a 
U.U4 


1 


2(a). pt(ji) > 125 GeV 


89 


67 


55 


21 


41 


2(b). Pt (j 2 ) > 50 GeV 


94 


93 


92 


76 


84 


3. 2 fe-tags 


22 


0.4 


6 


9 


0.3 


4. fir > 225 GeV 


56 


< 2 


< 5 


< 3 


< 10 (95% c.l.) 


v cxp (ioo fb- 1 ) 


7 


< 2.4 


< 2.7 


< 8.8 


< 177 (95% c.l.) 



Table 5: Same as Table 4, with an alternative set of requirements including 2 6-tagged jets. 



statistics, with only 7 signal events expected in a 100 fb" 1 data sample. 

To summarize, our analysis indicates that, for the MSSM parameters at the benchmark 
point, the signature (24) of the split-stop spectrum can be discovered at the LHC. The chosen 
BP is typical of the golden region, and this conclusion should generally hold as the MSSM 
parameters are varied away from the BP, scanning this region. There are, however, several 
exceptional parts of the parameter space where the observability of this signature could be 
substantially degraded. These include: 

• Large m 2 region: The I2 production cross section drops rapidly with its mass, see 
Fig. 3, suppressing the signal rates; 

• Small 9 t region: While non-zero 9 t is required in the golden region, values as small as 
0t = 7r/15 are allowed (see Fig. 2). The branching ratio Br (£2 —> Zt\) is proportional 
to sin 2 26* t , see Eq. (21), and the event rate is suppressed at small 9 t ; 

• Small ii-LSP mass difference: The absence of hard jets in this case would make the 
signal/background discrimination more difficult. 

In these special regions, observing the signature (24) may not be feasible at the LHC. 
These limitations should be kept in mind when theoretical interpretation of a search for the 
signature (24) is given. 

5 Alternative Interpretations of the Z + 2jh + fir + X 
Signature 

Unfortunately, observing an excess of events in the channel (24) at the LHC does not prove 
that the decay t 2 — > t\Z is occuring. Even within the MSSM, this is not the only possible 
interpretation of such an excess. The simplest alternative interpretation is stop or sbottom 



16 



t/b 



h 

(a) 



X 



o/± 
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Figure 5: Cascade decays in the MSSM leading to the Z + 2j + fir signature: (a) the chain 
characteristic of the golden region; (b) an alternative chain. 



production, followed by a cascade decay containing a b quark and a Z boson from a neutralino 
or chargino decay: x°j ~^ Xi% (due to the higgsino components of the neutralinos) or xf ~^ 
Xi Z. Can this alternative interpretation be ruled out based on data? 

One useful input for discriminating between these two interpretations is whether a signal 
is observed in a search identical to the one presented in Section 4, but requiring that the jets 
not be 6-tagged. If the signal is due to t 2 — > t\Z, all signal events contain energetic b quarks, 
and the number of events in this search would be zero if b tagging were perfect. The actual 
number of expected events under realistic conditions can be deduced from the error rate in b 
tagging, which can be measured elsewhere. If the signal is due to x°j ~^ Xi% or X2 ~ > Xi 
it does not have to be associated preferentially with third-generation squark production, and 
the number of events without b tags could be substantially larger than this expectation. This 
argument could be used to rule out the t 2 -^t\Z interpretation. Unfortunately, however, it 
cannot be used to confirm it: the pattern consistent with the t 2 — > t\Z interpretation may 
also appear if the events are actually due to chargino or neutralino decays, provided that 
the first two generations of squarks are substantially heavier than their counterparts of the 
third generation and their production cross section is suppressed. A direct measurement of 
the squark masses could break this degeneracy. If the first two generations of squarks were 
found to be light, but no signal is seen in Z + 2j + fir with non-6 jets, the "split stop" 
interpretation of the signal (24) would be preferred. 

A more direct way to discriminate between the two interpretations would be to study the 
distribution of the events as a function of the Z-jet invariant mass Sjz = {pj +Pz) 2 - This 
strategy is the same as the recently proposed method of discriminating between SUSY and 
alternative theories with same-spin "superpartners" [32, 33], but in this case it is applied 
to distinguishing two processes within the MSSM. Consider the Feynman diagrams corre- 
sponding to the two interpretations of the signal, shown in Figure 5. In the case of t 2 —*t\Z 
decays, the Z and the jet are separated by a scalar (stop) line, and their directions are 
uncorrelated. In the case of chargino or neutralino decays, the Z and the jet are separated 
by a fermion line, and spin correlations between their directions are possible. Unfortunately, 
in the neutralino case, no such correlations occur, because of the non-chiral nature of the 
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XiX°jZ coupling [33]. In the chargino case, however, the coupling has the form [1] 

CS(l-75) + Cg(l + 7 5)]fc> (25) 



2c 



w 



where 



C$ = UrjJv + -U* 2 U J2 . (26) 

Here, U and V are the rotations of the negatively-charged and positively-charged charginos, 
respectively, required to diagonalize the chargino mass matrix. Since in general U ^ V, the 
couplings (25) are generically chiral. The stop-bottom-chargino coupling is also generically 
chiral: it has the form 

b{A%P R + A%P L )x c ^ (27) 
which can be equivalently rewritten as 

-UA%P R + A%P L )b%. (28) 

Here 

S/3 

A% = ^U i2 R] 1: (29) 

C/3 

where R l is the matrix diagonalizing the stop masses: {ti,t2) T = i?'(ti,t^) T . Squaring the 
matrix element M. for the decay tj — > b + Z + xt , see Fig. 5 (b), and summing over the 
final-state polarizations yields 

£ l-^l 2 oc (\Af 2 \ 2 - \Af 2 \ 2 ) (|Cf 2 | 2 - |Cf 2 | 2 ) s bZ + const, (30) 

pol 

where the constant terms do not depend on Sbz, narrow- width approximation for xf nas been 
used, and the b quark mass was neglected. The charge-conjugate decay t* ^b + Z + xI has 
the same asymmetry. Observing a linear dependence of the event rate on s^z would provide 
clear evidence against the interpretation of the signal in terms of the process in Fig. 5 (a). Of 
course, in a real experiment, the asymmetry would be partially washed out by combinatoric 
backgrounds, as well as possible non-chiral decay chains containing the same final state. A 
detailed analysis of the observability of the correlation in Eq. (30) is beyond the scope of 
this paper. 

While our analysis so far focused on the decay t 2 ^ fiZ as a signature of the MSSM 
golden region, there are two other, closely related decays that are also characteristic of this 
region: 

t 2 ^b L W + , b L ^hW-. (31) 
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For example, at the benchmark point used for the analysis in Section 4, these two decays have 
branching ratios of 15% and 43%, respectively. Stop or sbottom pair-production followed by 
these decays leads to a signature 

W + 2 3b + fk + X. (32) 

This signature is complementary to the Z + 2jb + ftr + X signature studied above. On the 
one hand, it suffers from higher backgrounds, since the W cannot be fully reconstructed 
in purely leptonic channels. On the other hand, its interpretation within the MSSM is 
somewhat cleaner. The leading alternative interpretation of the signature (32) is that the 
W's are produced in chargino — > neutralino decays. But the chargino-neutralino coupling 
is chiral, and the directions of the W and the associated jet are correlated. If the W is 
sufficiently boosted, this will result in an observable linear dependence of the cross section 
on sij = (pe +Pj) 2 , where £ is the lepton daughter of the W [33]. If, on the other hand, the 
W is produced in decays of scalars, such as the processes (31), the distribution of events in 
S£j should be flat. 

To summarize, even if the MSSM is assumed to be the underlying model, the interpreta- 
tion of events with vector bosons associated with jets and missing E T is not unambiguous. 
Careful comparisons of the rates with and without b jets, as well as the distribution of events 
in vector boson-jet invariant masses, would be required to remove the ambiguity. This may 
take considerably more data than the discovery of an excess over the SM backgrounds in 
these channels. 

If the MSSM is not assumed from the beginning, the question of interpretation becomes 
even more confusing. For example, in the Littlest Higgs model with T-parity [34], Z bosons 
can be produced in the decay W\ — > ZB H , due to the mixing between the 577(2) and U(l) 
heavy gauge bosons. A similar decay involving the Kaluza-Klein states of the SU{2) x U(l) 
gauge bosons can occur in models with universal extra dimensions (UED) [35]. Again, 
a careful study of spin correlations would be necessary to disentangle these possibilities. 
Understanding the nature of such correlations in various models is an interesting direction 
for future work. 

6 Conclusions 

In this paper, we discussed an LHC signature of the MSSM characteristic of the "golden 
region" in the model parameter space. The advantage of this signature is that it directly 
probes the features of the stop spectrum that are dictated by naturalness and the Higgs 
mass bound. Experimentally, the signature is not straightforward, but the results of our 
simulations indicate that it should be within reach at the LHC. 

Given the strong theoretical motivation for the signature discussed here, we encourage ex- 
perimental collaborations to perform a more detailed study of its observability. The analysis 
of this paper relied on a set of simple rectangular cuts, and no systematic procedure to opti- 
mize the cuts was employed. It is very likely that a better algorithm for signal/background 
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discrimination, perhaps using modern data analysis tools such as neural networks or decision 
trees, would significantly enhance the reach. On the other hand, it should be noted that we 
ignored systematic uncertainties on the background rates in our reach estimates, and that 
no fully realistic detector simulation was attempted. 

If the first round of the LHC results points towards an MSSM-like theory, obtaining 
experimental information about the stop spectrum, and in particular testing whether the 
"golden region MSSM" hypothesis is correct, will become an important priority for the LHC 
experiments. An indirect way to shed some light on this issue by identifying the stop loop 
contributions to the Higgs production cross section and mass has been recently proposed by 
Dermisek and Low [36]. This is complementary to the direct probe explored in this paper. 
It would be interesting to explore other experimental consequences of the golden region 
hypothesis. 
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